Analysis of a microscopic stochastic model of microtubule dynamic instability 
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A novel theoretical model of dynamic instability of a system of linear (ID) microtubules (MTs) in 
a bounded domain is introduced for studying the role of a cell edge in vivo and analyzing the effect 
of competition for a limited amount of tubulin. The model differs from earlier models in that the 
evolution of MTs is based on the rates of single unit (e.g., a heterodimer per protofilament) trans- 
formations, in contrast to postulating effective rates/frequencies of larger-scale changes, extracted, 
e.g., from the length history plots of MTs. Spontaneous GTP hydrolysis with finite rate after poly- 
merization is assumed, and theoretical estimates of an effective catastrophe frequency as well as 
other parameters characterizing MT length distributions and cap size are derived. We implement 
a simple cap model which does not include vectorial hydrolysis. We demonstrate that our theoret- 
ical predictions, such as steady state concentration of free tubulin, and parameters of MT length 
distributions, are in agreement with the numerical simulations. The present model establishes a 
quantitative link between microscopic parameters governing the dynamics of MTs and macroscopic 
characteristics of MTs in a closed system. Lastly, we use a computational Monte Carlo model to 
provide an explanation for non-exponential MT length distributions observed in experiments. In 
particular, we show that appearance of such non-exponential distributions in the experiments can 
occur because the true steady state has not been reached, and/or due to the presence of a cell edge. 

PACS numbers: 87.16.Ka, 82.35.-x, 05.40.-a 



I. INTRODUCTION 

Microtubules (MTs) are intracellular polymers which 
provide a part of the cytoskeleton and are responsible for 
many cell functions including division, organelle move- 
ment, and intracellular transport. A cell is a living ob- 
ject, and as such it has to constantly adjust to and com- 
municate with a changing environment. For this purpose, 
MTs possess a property called dynamic instability, which 
enables them to promptlv switch between two modes, 
growth and shortening [ll, S Q ■ This is achieved through 
MT having a stabilizing cap which keeps the MT from 
disassembling. The MT tends to depolymerize when the 
cap is lost Q,y,|j,|^. The cap gradually hydrolyzes and 
becomes unstable as well, and so for the MT to survive 
it has to grow to renew its cap. 

The existence of a GTP cap at the end of MTs |j| and 
the phenomenon of dynamic instability 01 were discov- 
ered in the early 1980s. Hill and Chen used a Monte Carlo 
approach to simulate this behavior J2, |y|, employing a 
representation of an MT in which its cap could consist of 
many units (heterodimers). One of the main outcomes 
of their work was a suggestion that a two-phase (cap, 
no cap) model of dynamic instability, based only on ob- 
servable macroscopic rates of phase and length changes, 
was sufficient to understand the behavior of the ensem- 
ble of MTs (cf. 0, figs. 4-6]). This approach has been 
prevalent since then in modeling the behavior of an en- 
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semble of MTs H H i, ES El E El El • One 

exception was the theoretical model of Flyvbjerg et al. 
E^E3 which considered an elegant theoretical model of 
the GTP cap dynamics. This model was based on micro- 
scopic constitutuve processes of spontaneous and vecto- 
rial hydrolyses inside the MT, and fiuctuating growth of 
the cap size. Quite recently many detailed models of a 
single MT began to emerge El El El IH US which try 
to incorporate biological details observed recently due to 
advances in the experimental techniques. In particular, it 
is now known from the experiments that the tips of MTs 
can have geometrical configurations typical to growing 
and shortening MTs, which differ from one another (e.g., 
01). This is closely related to the idea of the structural, 
and not necessarily a GTP cap [131 , when due to ten- 
sile stresses inside the elastic body of an MT, its shape 
deforms from a cylinder towards the tip. 

In both in vitro and in vivo experiments, dynamics of 
MTs have been observed under a large variety of physical 
conditions and in various chemical environments. Much 
data has been accumulated including parameter values 
describing the MT dynamics and length distributions. 
Can these values be predicted based on the conditions of 
the experiments? How would change in ambient condi- 
tions or the presence of spatial constraints affect observ- 
ables? These questions are difficult, if not impossible, 
to answer using the models with postulated observable 
(macroscopic) rates. 

In this paper we analyze a novel model of MT dynam- 
ics in a finite domain bounded by the cell edge, which 
involves competition of individual MTs for tubuHn. The 
model is based on a linear ID approximation of a MT 



structure. We consider the role of the boundary and 
extend the model to incorporate finite hydrolysis. Our 
model is difi^erent from earlier works [TJ ll^ addressing 
the role of the edge in that we explicitly consider the 
concentration dependence of the dynamic instability pa- 
rameters, as well as a competition for a limited tubulin 
pool. 

Namely, we use a generaHzation of a microscopic model 
of MTs introduced in J23| ■ Instead of postulating macro- 
scopic rates J2; HI or deducing them from numerical sim- 
ulations [y|, we estimate them theoretically from basic 
microscopic rates of (de) polymerization and hydrolysis 
of a single unit (which can, but does not have to be in- 
terpreted as a heterodimer). This results in a higher- 
resolution theoretical model which may be more suitable 
for today's higher-resolution experiments, and can par- 
tially address the above questions. 
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One of the main results of the paper is an estab- 
lishment of a link, by using analytical formulas, be- 
tween microscopic parameters, describing polymeriza- 
tion/depolymerization and hydrolysis of individual units, 
and macroscopic (observable) characteristics of the MT 
dynamics and ensemble. We demonstrate how to ap- 
proximate macroscopic steady state behavior of MTs us- 
ing microscopic rates and vice versa, extract microscopic 
rates from macroscopic behavior. Hence, it makes it pos- 
sible to analytically and quantitatively predict the effect 
of changes in microscopic parameters on observable fea- 
tures, as well as to deduce microscopic changes from ob- 
served changes in macroscopic behavior, when relevant 
geometry and chemistry is taken into account. 

We also use a computational Monte Carlo model to 
provide an explanation for non-exponential MT length 
distributions observed in experiments. In particular, we 
show that appearance of such non-exponential distri- 
butions in the experiments can occur because the true 
steady state has not been reached, and/or due to the 
presence of a cell edge. 

The paper is organized as follows. The conceptual 
model and its computational implementation are pre- 
sented in Sectionin Next, in Section lTTTl we develop a cap 
model where the cap can have any number of (sub)units, 
which can be single heterodimers - see Figure Qp. This 
cap model differs from previous models in that it does 
not involve vectorial or induced hydrolysis. Using this 
model we derive approximate expressions for observable 
rates. We then describe in Sections ll VI and FVl a quantita- 
tive theoretical analysis of a lower-resolution model with 
the cap being treated as a single unit - see Figure Q]^. 
The influence of the cell edge is also studied there. Sec- 
tion describes balance between polymerized and free 
tubulin in a bounded domain with a fixed total amount 
of tubulin present. Finally, we discuss and summarize 
our findings in Sections I VII and fVIII 



Figure 1: Schematic representation of the model. 



II. MODEL DESCRIPTION, PARAMETERS 
AND NOTATIONS 

In this section we describe a basic model of dynamics of 
MTs. We consider a domain of size L^x LyX L^ with Nn 
available nucleation sites for MTs in its center. Nn is the 
maximal number of MTs - cf. Table |l| For simplicity, in 
our study of the role of the boundary (e.g., cell edge) we 
assume that all MTs have an identical maximal allowed 
length. All MTs grow from nucleation sites (there is no 
spontaneous nucleation) , and MTs grow at one end (usu- 
ally the so-called plus end) only. There is a fixed amount 
of total tubulin in the domain. This tubulin is present in 
two forms: free tubulin in the solution and polymerized 
tubulin constituting the MTs. Free tubulin is taken up 
by growing (polymerizing) MTs and is released back into 
the solution by shortening MTs. In general, free tubuHn 
(Tu) diffuses inside the domain. In this paper we assume 
that the diffusion of free Tu is fast and does not lead to 
a diffusion-limited reaction rates. This is in agreement 
with J2J|. Moreover, we assume uniform concentration 
of free tubulin throughout the domain which implies in- 
stantaneous diffusion^ (For the studies of the effects of 

tubulin diffusion see (23,l2a|-) 

In this paper a MT is represented at each moment in 
time in the form of a ID straight line consisting of a cer- 
tain number of units of a predefined length. Each unit 
belonging to a MT can be in either a growth-prone state 
or a shortening-prone state. We will refer to them as 
GTP (or T) state or GDP (or D) state respectively. All 
free tubulin is assumed to be in a T state. When a unit 
joins the MT it is initially in a T state. The probabil- 
ity that the internal units have hydrolyzed (transformed 
to a D state) increases with time. When MT disassem- 
bles (shortens) these D units, upon becoming terminal, 
have higher probability to disassemble and return to the 



Symbol 


Definition 


Dimensions 


c 


concentration of free tubulin 


pM 




critical concentration of free tubulin 


fiM 


Ctot 


total concentration of tubulin 


fiM 


kgT 


pseudo-first order rate of adding a unit on top of a terminal T unit 


/^M-ifi-i 


kgD 


pseudo-first order rate of adding a unit on top of a terminal D unit 


fiM'^s-'' 


X 


parameter of exponential distribution (number of units) 


- 


£ 


characteristic cap size (number of units) 


- 


m 


mean length (number of units) of MTs 


- 


n 


coarsened step size in the cap model (number of units) 


- 


P 


(n-l)/(^-l) 


- 


K, 


rate of the edge-induced catastrophe 


S-' 


Kh 


rate of hydrolysis (transformation to D state) of internal T units 


S-' 


KgT 


rate of adding a unit on top of a terminal T unit 


S-' 




effective rate of growing by one unit in growth phase 


S-' 


J^qT 


rate of growing by n units in growth phase 


s-^ 


KgD 


rate of adding a unit on top of a terminal D unit 


S-' 


K„ 


nucleation rate of a MT 


S-' 


KsT 


rate of depolymerization of terminal T unit 


S-' 


T^eff 


effective rate of shortening by one unit in growth phase 


s-^ 


T^obs 
"■sT 


rate of shortening by n units in growth phase, = catastrophe frequency 


s-1 


KsD 


rate of depolymerization of terminal D unit 


S-' 


L 


maximal length (number of units) of MTs in domain with upper bound 


- 


J^x^ -^y, J^z 


domain sizes in the numerical simulations 


m 


Mg{l) 


number of MTs of length I (number of units) in growth phase 


- 


Ms{l) 


number of MTs of length I (number of units) in shortening phase 


- 


No 


number of free nucleation seeds 


- 


Nmt 


number of MTs, < Nn 


- 


N„ 


total number of nucleation seeds 


- 



Table I: Notation highlights. The dash (-) means that the considered parameter/variable is dimensionless. 



solution than the terminal T units. Upon return to the 
solution they immediately switch to the T state. The ter- 
minal T unit does not hydrolyze but can with a certain 
probability depolymerize (drop from MT end). Incorpo- 
ration of a new unit at the MT tip triggers the hydrolysis 
process of the previously terminal unit. This description 
seems appropriate in view of J3|. 

The dynamics of the MTs is determined by five mi- 
croscopic rates, KgT (Kgo) and Kgx (Ksd), which are 
the rates of MT growth and shortening (i.e. adding one 
more unit from the solution on top of the current ter- 
minal unit, and losing this current terminal unit to the 
solution) when the terminal unit is in state T (D), and 
the hydrolysis rate Kh of the internal units which are in 
state T. If the terminal unit has to hydrolyze in order to 
depolymerize (and its hydrolysis rate is not faster than 
that for the internal units) then Kgx < Kh- 

For numerical simulations, shortening rates are taken 
to be independent of c while growth rates are assumed 



proportional to c at the location of MT tip: 



KgT,gD = kgT,gDC. 



(1) 



Such specific dependence of the growth and shortening 
rates on c, though, is not required for many of the theo- 
retical results we report. 

When the MT reaches boundary of the domain it is 
not allowed to grow any more and will eventually lose its 
terminal unit initiating with certain probability a short- 
ening phase. There are two more rates at the domain 
boundary of importance in the model: rate Kn of nucle- 
ation from existing seeds and a rate K^ of edge-induced 
catastrophe, which can also depend on c. Appendix El 
contains a brief description of a numerical algorithm we 
used in our simulations. 

In what follows, we will impose restriction on a max- 
imal length of a MT (upper bound, e.g., due to a cell 
edge). We will call zero a lower bound. 




30 40 50 
time (s) 



Figure 2: Length history plots of three arbitrarily selected 
MTs. Zero time here corresponds to 80 seconds from the 
start of the simulation. Here Kh — 10, ctot ~ 10 and other 
parameters are listed in Table |ll| See also Fig. |^ 



A. Observables 

The standard experimental observables describing the 
dynamic behavior of a single MT are derived from the 
MT length vs. time plot. Typical length history plots of 
MTs are shown in Fig. [2l Due to the two-state nature of 
the tubulin units inside the MT, the fluctuations in this 
length may be large and even in the (macroscopically) 
steady state each MT can repeatedly change its length 
all the way from zero to some characteristic length, or to 
the boundary. If no boundary is present and if free tubu- 
lin concentration stays high enough (if Knr and Kgjj are 
high enough) MTs can grow unbounded [3,l2j,l23- From 
a sawtooth-like evolution of an MT length four param- 
eters can be extracted: the velocity/rate of growth, the 
velocity/rate of shortening, the average time of growth 
before switching to shortening, and the average time of 
shortening before switching to growth. The inverses of 
these times define the so-called catastrophe frequency 
and the rescue frequency, respectively. Note that brief 
growth or shortening intervals may pass unnoticed in the 
analysis of experimental data. Some models of dynamics 
of MTs use these four parameters as g iven constants for 
constructing analytical solutions PiUj ignoring their mi- 
croscopic origin. In Section ITTTI we use microscopic rates 
to derive the observable growth velocity and catastrophe 
frequency instead of setting them from the beginning. 



III. CAP MODEL 

Carlier and colleagues [J, |23, |2^ have provided exper- 
imental evidence that the GTP cap of the MTs is not 
restricted to the units at the very tip. This suggests that 
the hydrolysis is not instantaneous, a conclusion also sup- 



ported by work on yeast tubulin |3(l. l3l] | . To incorporate 
this feature into our approach we develop a new model for 
the cap. In this section we show that using this model, 
based on the underlying microscopic laws, one can pre- 
dict the observables: the catastrophe frequency and the 
velocity (rate) in the growth phase of MTs. In other 
words, the microscopic laws of the MT dynamics, gov- 
erning single unit polymerization/depolymerization and 
hydrolysis, can be related to macro-scale (observable) dy- 
namics of an MT. 

When Kh < oo, the MT cap in our model consists 
mainly of T units (cf. Fig. Qp) and possibly of a few 
D units and has some characteristic length (number of 
units) £ > 1. When Kh — » oo then £ = 1 because only 
the terminal unit is not allowed to hydrolyze and it is 
in a T state. Our approach is based on coarsening the 
resolution in the growth phase so that only blocks of the 
order of cap size, n ^ £, are added. By catastrophe we 
understand the loss of the cap. In what follows we will 
establish a connection between coarsened "observed" rate 
constants K°^gj, of growing or shortening in the growth 
phase, and the original rates KsT,gT and K^. There is no 
need to rescale KsD,gD, if the free tubulin concentration 
is not too low. Therefore, we only consider a model for 
the MT cap and don't alter the rest. 

Let us consider the cap model in detail. In what follows 
we neglect the fluctuations in £ due to randomness in hy- 
drolysis, and we assume that each T unit is hydrolyzed 
after staying an internal unit for a time At/i = 1/Kh- 
After the rescue or nucleation event occurs the cap be- 
gins to grow. It has time Aig + At/j to elongate, where 
Aig = l/{KgT + Kst) is the time of a single-unit step 
in the growth phase. After that its average length re- 
mains constant (under assumption that no catastrophe 
occurs during this time). For the catastrophe to occur 
the cap should be lost, due to fluctuations in cap size 
and in the growth velocity [3,IMI^- In our analysis we 
consider two scenarios for cap loss: (i) roughly half of 
the cap is lost due to random nature of MT's growth, 
and the second half gets hydrolyzed during this time, or 
(ii) the whole cap is lost due to random fluctuations in 
MT growth. Keeping in mind that the terminal unit can- 
not be lost as a result of hydrolysis, in description (i) we 
require that {£ + l)/2 and {£ — l)/2 units are lost due to 
fluctuations in MT growth and propagation of hydrolysis 
front, respectively. We deflne 

n ~ \ 



P = 



i-\ 



(2) 



In case (i) n ~ (£ -f l)/2 and p = 1/2, while in case 
(ii) n = £ and p = 1. It is important to stress that p 
is introduced as a flxed parameter set a priori, based on 
the scenarios of cap loss similar to (i) and (ii). The two 
descriptions (i) and (ii) determine the duration of the 
coarsened step 



Atf' = — . 



1 



gT 



1 1 
r- = pAth+Ata = p \ 

K°^' ^ '^ ' ' Kh KgT 



KsT 

(3) 



For a given n, in order to rescale/coarsen the dynamics 
in the growth phase we require that both the average 
velocity and the diffusion coefficient of the MT tip remain 
unchanged. For a random walk on a line, with probability 
p to jump to the right and q = 1 —p to jump to the left, 
the average velocity is w = {p—q)Ax/At and the diffusion 
coefficient is D = 2pqAx^/At. Here Ax is the step length 
and At is the time per step. In case of the original walk 
p = Pg ~ KgrAtg, Ax = 1 aud At ~ Atg, while in 
case of the rescaled walk p = K°i^At°'"', Ax = n and 
At = Atg^^ is given by eq. IJSJ. After introducing the 
effective rates of adding or losing one unit (as opposed to 
one block), 



and using conservation of v and D we obtain that 



KfJ - K%f 



KfjK^'n 



K, 



gT 



K 



sT 



KgTK.T 



-eff 



^eff 



Kt+Kt "^aT 



K, 



K 



sT 



(4) 



(5) 



(6) 



From eqs. (EJ and l|l|| n/{K''JJ + K^JJ) can be found 
and substituted into eq. lO, resulting in 



KrKi' = P, 



where 



(i = 



K„tKst 



l + p- 



K, 



gT 



(7) 



(8) 



where the term K'^Jf Ath approximates the number of 
added units after the beginning of a growth phase, before 
the hydrolysis front starts moving. In fact, eq. lfT2|l can 
be used as a definition of £ and n can be found from eq. 
(j2Il. Then there is no need for eq. JSJ as eqs. J^l, © 
and lfT2)l form a closed set of equations. This approach, 
however, leads to a cubic equation for K'^U , while in the 
above approach we need to solve a quadratic equation, 
which is much simpler. Nevertheless, the definition of 
I through eq. ltT2)l seems to work better in the limit 
of KgT -^ 0. Namely, substituting eq. lfT2)l into eq. 
(|2j) yields n = 1 + pK'^fp' /Kh and the only non-negative 
solution of this equation together with eqs. Q and I©, 



in the limit KgT 



0, is KlfJ 



0, K%f = KsT and 
n = 1. Indeed, it seems reasonable to postulate that 
n ^ 1, i.e., there is no rescahng, when KgT -^ 0. This 
does not follow from eqs. (jHJ and ltTT|l . which lead to 
n ^ 1 + pKst/ Kh > 1 instead. 

Using the above developments, it is possible to derive 
scaling behaviors of various quantities as functions of, 
e.g., c and K^. For example, substituting eqs. Q and 
II12|I into eq. (jHJ and using eq. Q together with eq. Q 
leads to 






1 



KgTK, 



Ky 



KfJ + {n 


-l){KgT+KsT) 


KsT 
7l2 ' 


KgT > KsT 


KsT 


1 



(13) 



After solving Eqs. ^ and Q and choosing only positive 
solutions we obtain that 



Kl 



-eff 


KgT ~ KsT + \/{KgT ' 


- KstY + 4/3 


gT 


2 




ff _ 


-KgT + KsT + ViKgT- 


- KstY + 4/3 



(9) 



(10) 



Expression (O displays some expected features. Namely, 



when Kh 



K 







K 



•ff 
gT 



eff 

gT 



KgT and when A'/,, 
{KgT - KsT + \KgT - Kst\)/2 



KgT - KsT, KgT > KsT 

0, KgT < KsT 

-eff 



. In general, max(0, KgT 



Kst) < K^^' < KgT. 

Eq. ^ combined with eqs. lO and (llOjl yields 



KgTKsT ^{KgT-KsT)^+^l3 



K. 



gT 



K 



sT 



P 



(11) 



and eq. ^ can now be used to determine L On the 
other hand, i can be approximated as follows: 



Klf/Ath + 1 



K 



eff 

gT 



(12) 



so that the catastrophe rate (frequency) K°!f scales as 
n~2 ex £~2. This is characteristic of diffusive scaling 
because time to catastrophe is determined by diffusive 
movement of the hydrolysis front relative to the MT tip. 
If free tubulin concentration c is not too small, then eqs. 
itTTjI . © and 1^1 yield n ex c and hence K°^^ ex c"^, 
which is in at least qualitative agreement with previous 
predictions [isl^Sl- The scaling Kfj? ex n"^ might have 
been postulated, as well, which would have led us to an 
additional version of the solution of the cap model. 

To provide another scaling example, let us now assume 
that Kh is small and KgT ~ KsT- Then from eq. (|Hl it 
follows that /3 « KhKsT/{2p) oc Kh. If /3 > {KgT - 
KstY , i.e., if Kh is not too small, then eq. Hll|l yields 
n ex Ky^ ' and hence I ex K^^ ' as well. Notice, that 
it is the same scaHng as derived for actin polymers (34 
eq.3]. 



IV. ENSEMBLE DYNAMICS OF 
MICROTUBULES 

In this section we treat cap as a single effective unit - cf. 
Fig. m Thus the model essentially reduces to the two- 
phase model proposed in ja|. First, we rederive length 



distribution of MTs, which are known in the Uterature. 
In particular, we study the role of upper bound (e.g., 
cell edge) . We use these results for analyzing in the next 
section competition for a finite tubulin pool. We also 
consider the steady state critical concentration of free 
Tu. 

In what follows, we use either the discrete or the con- 
tinuous description of MT dynamics, whichever is con- 
venient. We assume that the continuous model provides 
a good approximation of the discrete model. The con- 
tinuous approach was discussed in [3, |^ ISa| while an 
analogous discrete approach was developed in 0, uM- 
Following [3 we write down the equations for length dis- 
tributions of MTs in growth and shortening phases in the 
form 



dtMa 



^KrMg 



KgoMs 



-JJ^ 



KlUdiNL 



(14) 



consider a steady state. Adding up these two equations 
then leads to di{—K'gip'Mg + KsdMs) = 0, which means 
that the spatial derivative of the flux (of the MT tips, 
considered as random walkers) is zero meaning that the 
flux is uniform. However, in the closed system this flux 
must be zero and hence 



M, 



K 



gT 



K. 



M„ 



D 



(19) 



It follows that eqs. lO, ^^ and (QHl still hold inside the 
domain, except for A not being necessarily positive, which 
is in agreement with previous work [ll|. This means qual- 
itatively that there might be a steady state distribution 
of MTs in which most of them are close to the upper 
boundary (e.g., cell edge), while only a few are short. 



ftM, 



K°^'Mg - K, 



gD 



Ms + KsodiMs 



(15) 



where Mg^s{l,t) are densities of MTs of length / at time 
t, in the growing (g) and shortening (s) phases. 

Equations Ijl4|l and l|15|l can be used to describe reg- 
ular diffusion with drift, if we don't distinguish between 
the phases (see also [1^1) • However, it is important to 
stress that these equations don't have diffusion terms for 
Mg^s{l,t) and hence switching phases back and forth is 
the only mechanism of spreading of these distributions 
present. This is in agreement with our simulations in the 
case of instantaneous hydrolysis of internal units. Notice 
that diffusion terms are used in Q- 

First, consider a semi-infinite domain. Equations ifnji 
and itlfijl with the boundary condition Mg^s{l = 00) = 
have the following steady state solution 



M„ 



Ae-'/\ 



where 



M. 



1 
A 



K 



eff 



KsD 



rrobs 

J^eff 
^^gT 



K 



gD 



K. 



sD 



(16) 



(17) 



(18) 



The necessary condition for the existence of a steady 
state in the case without a boundary is given by A > 0. 
The prefactor A is a normalization coefficient which 
depends on the total number of MT nucleation seeds 
present and on the nucleation probability. 

We now add a constraint limiting the maximal length 
of MTs to be L. MTs cannot become longer due to a 
barrier, for example a cell edge, as is often the case in 
vivo, especially when the cell is in the interphase and 
the MTs are relatively long. For simplicity we assume 
that L is identical for all MTs. We still can use eqs. 
II14|I and II15|I inside the domain, for Q < I < L, and 



A. Critical concentration of free tubulin 

Let us consider the Hmiting case 1/A = (cf. eq. CHI)- 
This defines the upper limit of the concentration of free 
tubulin c|^ at which the steady state in the semi-infinite 
domain still exists (cf. PI)- Let us use eq. Q and define 



^sT 



IKstK, 



TJ^sD 



'igT 



kgrkgD 



(20) 



Because the MTs in the growth phase are less likely to 
shorten than are MTs in the shortening phase, a < b. In 



general. 



-eg 



G \a. 



b] . Notice that slowdown of hydrolysis 



reduces c'^{Kh). When hydrolysis is instantaneous then 
c reaches its maximal value c^(oo) = b. When there 
is no hydrolysis at all then c'^{Kh) reaches its minimal 
value c^(0) = a meaning that the average growth rate 
in this case, kgTc'^{0) — Kst, is zero. 

One can get a scaHng estimate of c^ if it is far enough 
from both a and b. Assume that KgT ^ KsT and 
KgT > Kh. Using eq. (Pi, then K^^,^ w KgT oc c. Now, 
p is of order of 1, so that from eq. (jHJ and eq. lfTT|l it fol- 
lows that /3 ~ KstKii <C K^j, and n ~ Kgx/Kh oc c/ Kh, 
respectively. Hence Kfj? k. KsT/n^ oc (Khjcf ■ Substi- 
tuting these scaling relations into 1/A = and using eq. 
(jIHl yields c^ oc ^/K^. 

Notice that if there are no rescues, kgo = 0, then c^ 
is infinite (see also [23|) and unbounded growth can not 
happen. This is so because without rescues the MT de- 
polymerizes completely after the catastrophe, no matter 
how long it was before. 



V. COMPETITION FOR TUBULIN 

In Section ITvl we have shown the existence of a steady 
state distribution of MTs inside a domain. It is conceiv- 
able that by sensing and controlling free tubuHn concen- 
tration and the number of MTs the cell regulates MT 



dynamics, as suggested in |2M l3^. In what follows we 
show in detail how to determine a steady state concentra- 
tion of free tubulin c, which is the key to finding steady 
state characteristics of MTs in a closed system. The main 
goal of this section is to derive expressions for the average 
number of units per MT m, and a number of MTs Nmt 
as functions of c and the other parameters. These func- 
tions are needed to determine c from the conservation of 
total tubulin. Because m can not extend beyond the do- 
mains' boundary, and because Nmt is less or equal than 
a number of nucleation sites 7V„, the number of polymer- 
ized units in a bounded domain stays restricted as the 
amount of total Tu grows. 

In what follows we assume that the total amount of 
tubulin is constant and we consider bounded domain of 
volume V . We also assume instantaneous diffusion so 
that the concentration of free tubulin is uniform through- 
out the domain. Hence, 



where /f„ is a nucleation rate, which in general depends 
on c. The left-hand side of eq. ll2ljl describes the rate of 
production of available nucleation sites by completely de- 
polymerizing MTs. The first term represents those MTs 
which experience a catastrophe, while the second term 
stands for those MTs which are already in the shortening 
phase. Using iflfijl . ((TTIl . and assuming that A ;» 1, yields 



A( A',T + Kl^' )-KnNo. (25) 

In addition, approximating summation by integration as 
follows 



NMT=J2{Mg{l)+Ms{l)) 



1=1 



K 



-eff ' 



gT 



sD 



AX (26) 



results in 



Ntot = Nfree + mN MT , 



(21) 



where Ntot is a total number of tubulin units in the do- 
main and Nfree is a number of free tubulin units. By 
dividing this formula by V we obtain expression for con- 
centrations measured in micromolars (/iM) 



Ctot 



tjiNmt 
IQ-^NaV 



(22) 



where V is given in m^ and IfiM is equal to lO^^A^^ 
units per liter or 10~^Na units per m^; Na ~ 6.022 x 
lO^^moP is the Avogadro's constant. We will use this 
expression for studying MTs in cases of unbounded and 
bounded domains. 



Nmt 



Nn 



irobs 
^sT 



^gT 



XK,,{l + Kli/ /K,d) 



gT 



(27) 



Substituting eqs. JSjl, JS^I and (jTHIl into eq. (|22)|, and 
using dependence of the rates on concentration c, eq. 
CJ, yields an equation Ctot = F[c)^ at steady state. This 
equation relates free tubulin concentration c and total 
concentration Ctot, as a function of all given parameters. 



B. Bounded domain 



Unbounded domain 



It has been shown in Section |^ that the steady state 
distribution of MT lengths in this case is exponential as 
described by eqs. Ijl6|l and Ijl7|l with A representing mean 
length of a MT 



m ~ \. 



(23) 



To find Nmt for a given Nn, we use a balance equation 
for the number of available nucleation sites A^o = Nn — 

Nmt'- 



Kt¥MS 



l)+KsDMsil = l)=KnNo, 



(24) 



Here again we impose limitation on the maximal pos- 
sible length of MTs not to exceed L. In the case of a 
bounded domain a steady state solution of eqs. ifnji and 
II15|I can be calculated, and then m and Nmt are de- 
termined in a way similar to the previous case. Notice 
that the steady state does exist even if c > c^. Since 
< m < L and < Nmt < Nn always hold, the second 
term in eq. H22II is non-negative and bounded. Therefore, 
when c goes from to oo so does ctot ■ If the right-hand 
side of eq. lf22|l monotonically increases with c then there 
is a unique physically meaningful c for each ctot , at steady 
state. 

It is shown in Appendix El that now 



N 



MT 



1 + 



Kn 



Nn 



T^obs 
^sT 



K 



eff 

gT 



(28) 



(1 + K;^^/K,d){1 - e-VA)A -1- iKf^f/Ke)il 



KgD/KsD)e 



-L/\ 



and 



{1 + Kfj/K,D)X{X 



^eff 



e-^l\L + X))+L{Kl'r,J/K,){l+KgDlK,D)e 



-L/\ 



(1 + Klf//K,D)\{1 - e-^/^) + {KlfJ/K,){l + KgD/K,D)e-^/^ 

I 



(29) 
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Figure 3: Frequency of catastrophe K"!^" as a function of MT 
growth velocity Vg. Dots represent experimental data [lfl . l37l |. 
Solid and dashed lines correspond to p = 1/2 and p = 1, 
respectively, and were obtained using the model f Section ITTTt . 



Inset: Two upper curves on the right are K'^Jf jKgT and two 
lower curves on the right are K^tJ /K^j. . 



When A > and L ^ oo eqs. ^^ and ((22|| reduce to 
eqs. ll?7jl and (|^ . respectively. When K'^^^ — > and 
KgD -^ 0, then A i (so that A > 0), Nmt -* and 

-^ cx), then A t 
L, as expected. 



0. When K 



eff 

gT 



(so that A < 0), Nmt 



oo and Kgjj 

N,, and m -^ 



VI. DISCUSSION OF THE RESULTS 



A. Comparison with existing cap models 



The cap model presented in this paper differs from 
the approaches used in [2, IM lift 113 • First, we don't 
postulate catastrophe frequency and growth velocity, or 
derive them from numerical simulations, as was done in 
[3,0|- Instead, we analytically derive these macroscopic 
rates from small-scale rates (such as chemical rate con- 
stants). Second, we employ only spontaneous hydroly- 
sis and don't use induced, or vectorial, hydrolysis; both 
types of hydrolysis were used in (ift Il7l |. Our model 
agrees with the experimental data analyzed in [ifl. Il7l | as 
can be seen from the main panel in Fig. Specifically, 
the predicted dependence of the catastrophe frequency 
on MT growth velocity is in agreement with experimen- 
tal data. 

The cap dy namics in the model proposed by Flyvbjerg 
et al. (ift [13 is modeled by addition of tubulin from the 



solution to the MT tip. This addition (polymerization) 
is faster than the propagation of the induced hydrolysis 
front (low end of the cap). Therefore, the cap length 
would grow infinite were it not for the spontaneous hy- 
drolysis at some point inside the cap. When it occurs, 
the cap is redefined as an interval between this sponta- 
neous hydrolysis point and the MT tip. In this way, the 
average cap size can be kept constant at steady state. 
According to Flyvbjerg et al. [ift Il7l|. catastrophe oc- 
curs when this cap is lost, and it is postulated that the 
remaining GTP-Tu units below the cap are not capable 
of rescuing the MT. This assumption is made in order 
to allow for catastrophe to occur. Otherwise, in many 
cases the rescue would immediately follow the cap loss. 
While this picture is widely accepted, in our alternative 
model the picture is even simpler. We use only one type 
of hydrolysis and we don't need to make any additional 
assumptions. In our model, there is a hydrolysis front 
due to spontaneous hydrolysis of old enough units (see 
Section Unll . The velocity of this front is governed by the 
age of the units inside the MT and hence it is always 
approximately equal (with fiuctuations) to the growth 
velocity. Faster growth velocity leads to a longer cap, 
reducing the catastrophe frequency. 

Dilution experiments have shown that sharp reduction 
in the concentration of free Tu to low or zero values re- 
sults in collapse of the MTs after a certain delay. Impor- 
tantly, this delay is practically independent of the initial 
free Tu concentration J2a|. Flyvbjerg et al. explain this 
phenomenon by arguing that the dilution results in dom- 
ination of spontaneous hydrolysis which regulates the 
waiting time before the collapse. Therefore, this time 
is almost independent of the initial cap size. Our model 
yields the following simple explanation. When concen- 
tration of free Tu becomes very low, two events must 
occur for the cap to disappear. The terminal unit should 
be lost (rate Kst) and the next unit should hydrolyze 
(rate Kh). If this next to last unit is old enough to hy- 
drolyze then, with high probability, the rest of the cap 
has already hydrolyzed. 

Dilution experiments reported in [23| and cited in [Tfll 
determine the average waiting time before the catas- 
trophe as roughly 5-lOs. Therefore, in Fig. El we set 
Kh — KsT ~ 0.15s~^. If the shortening velocity (Ksd) is 
much larger than KsT , and if the loss of the terminal unit 
is conceptualized as a two-stage process of hydrolysis and 
then falling, the equality Kh = K^t would indicate that 
the hydrolysis rate of the terminal unit equals the hydrol- 
ysis rate of the internal unit. Notice that we were not 
able to fit the catastrophe frequency data if K^t and Kh 
were significantly different. In the dilution experiments 




Figure 4: Steady state concentration of free vs. total Tu 
in the unbounded (thin lines) and bounded (thick lines) 
domains. For a bounded domain maximal MT length L 
is roughly 1400 units. For a bounded domain c reaches 
its asymptotic value c^{Kh). Hydrolysis rate is K^ = 
10s~^. Full and dashed lines correspond to p = 1/2 and 
p = 1, respectively. Other parameters are Nn = 200, 
V = lO-^m^' = IO^VL, kgT = 5^M-^s-\ K,t = 5s-^ 



kgD = 0.5fiM-^s-\ KsD = 500s-\ Kr, 
and we use eq. Q. 



KgD, K, 



K, 



spatial resolution was about 0.25/iTO j^al; which is about 
30 heterodimers (per protofilament) , so that actual wait- 
ing time before losing the terminal unit (heterodimer) 
might be faster than the reported waiting time before 
the collapse. We used eq. QJ with kgT = 0.3^M~^s~^, 
however any value of kgT can be used, because it en- 
ters the formulas only through KgT = kgTC and there 
is no explicit c-dependence in the figure. The unit 
length is taken to be the length of one heterodimer of 
Tu, 8nm and hence lunit/s ~ 0.5 ^m/min. Therefore 
Vg {iim/min) « K'^JJ /2 {units /s). In the inset of Fig. 
we show the values of K^fp^ /K'^fJ and K^JJ jKgT- It is 
seen that for «„ > 0.2 one has K°X^ < Kllf < Kill and 



sT 



hence Vg ex K^^t ^ '^• 



^gT 



B. Competition for tubulin and the edge effect 

Here we study competition for a Hmited pool of free 
tubulin and combine it together with the cap model. At 
steady state, the dependence of free tubulin concentra- 
tion, c, on the total tubulin concentration, Ctot, is gov- 
erned by tubuHn mass conservation eq. lf22|l . The result- 
ing value of c, in turn, defines the dynamics of MTs and 
their ensemble distributions. 

In the case of an unbounded spatial domain, we have 
reproduced the prediction of the Oosawa and Kasai 
model [231 including existence of a critical concentration 
of free tubulin c^ (see thin lines in Fig. 31- The steady 
state cannot exist above this value. Depending on the 
choice of parameters, the transition from almost linear 



growth of c, c w ctot for low Ctot, to the asymptotic value 
c = c^ can be made sharp or smooth. (In Oosawa model, 
this transition is assumed to be sharp meaning that when 
Ctot < c^ there are no MTs and when Ctot > c^ all the 
excess tubulin above c^ goes into polymerized state.) 
When the probability of rescue is zero then c ^ — > oo, 
our model describes the situation considered in |33 |. 



If there is an upper bound on MT lengths, then at suf- 
ficiently high concentrations of total tubuHn, this bound 
inhibits further polymerization of MTs (edge effect). 
Hence, the steady state free tubulin concentration can 
rise above its critical value. This edge effect is demon- 
strated by thick lines in Fig. Eland was first discussed 
in J23|. Under our reference conditions, it is seen that 
for Ctot > 20^M the edge starts playing an important 
role. At sufficiently high ctot, the edge reestabHshes lin- 
ear growth of c with respect to Ctot- The implications 
of this effect are as follows. If Ctot is high enough, the 
MTs grow persistently up until hitting the edge which 
triggers their catastrophe. This is consistent with recent 
experimental observations of persistent growth of MTs 
in vivo J4Q|. From our model it follows that by changing 
the number of nucleation sites while keeping the total 
amount of Tu constant the cell could regulate the transi- 
tion between mitotic (short) and interphase (long) MTs. 

Next we compare in Table O Monte Carlo simula- 
tions (see Appendix EJ with the results obtained by 
using our continuous model. We choose large domain 
size so that MTs never reach the boundary and m ~ \ 
(cf. Eq. (l23|l ). We run simulations for different pa- 
rameter sets until steady state is reached and then de- 
termine free tubulin concentration, number of MTs and 
their mean length, and estimate the cap size. Recall 
that only MTs in growth phase have caps. Therefore, 
we estimate the cap size I of MTs in our simulations by 
l^ {1 + Kl{,^' /Ksd) ■ #T/Nmt, where #T is the number 
of polymerized Tu units in T state. We also use eq. (I19II . 
For the parameter values used in Table ^simplified for- 
mula for the estimated cap size £ ^ ^T/Nn has been 
used. 



TahlelTTIalso contains theoretical estimates correspond- 
ing to the simulated values discussed above. In addition, 
the table includes theoretical estimates of c?S and of the 
cap size based on eq. irL2|l . In most cases the simulated 
results He in between our two theoretical approximations, 
for p = 1/2 and p = 1 respectively, in agreement with 
the model description of cap evolution. These approxi- 
mations are given as two adjacent numbers in the cells of 
the table displaying our theoretical estimates. When, 
however, Kh becomes small, and c approaches a (eq. 
Ij20ll ). our approximations seem to consistently overesti- 
mate the number of MTs, Nmt- This should be improved 
by rescahng the rest of the rates, i^gu, i^sD, K^ and K^, 
which is outside of the scope of this paper. 
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Parameters 


Simulated results 


Theoretical estimates: n = [(. + l)/2, n = (. (see text) 


Kh 


Ctot 


iVn 


c 


m 




Nmt 




C 


A 


e. 




Nmt 


oo 


36 


2000 


31.32 


1416 


0.759 


1990 


31.62 


31.28 


1430 


1 


1 


1989 


oo 


28 


3000 


27.46 


117 


0.734 


2814 


-'■- 


27.48 


112 


1 


1 


2797 


oo 


28 


10'', T seed 


26.48 


92.7 


0.780 


9908 


-'■- 


26.53 


89.5 


1 


1 


9909 


10^ 


36 


2000 


31.36 


1401 


0.997 


1995 


31.60, 31.57 


31.25, 31.23 


1437, 1444 


1.0015 


1.0016 


1989 


100 


30 


2000 


20.78 


2788 


1.75 


1992 


21.05, 17.38 


20.92, 17.29 


2742, 3836 


2.00, 1.82 


2.03, 1.84 


1994, 1996 


30 


15 


2000 


12.7 


694 


2.69 


1969 


14.7, 11.5 


13.9, 11.2 


347, 1139 


3.16, 2.72 


3.22, 2.77 


1952, 1985 


10 


10 


2000 


7.51 


758 


3.94 


1974 


9.85, 7.48 


9.02, 7.18 


305, 857 


5.04, 4.12 


5.16, 4.20 


1945, 1980 


10 


5 


2000 


4.91 


35.1 


2.25 


1538 


-'■- 


4.95, 4.85 


21.9, 52.1 


3.06, 2.98 


3.20, 3.07 


1420, 1722 


3 


5 


2000 


3.94 


327 


5.66 


1954 


6.08, 4.59 


4.64, 3.99 


115, 312 


7.22, 6.09 


7.42, 6.21 


1868, 1952 


1 


4 


2000 


2.56 


443 


8.51 


1958 


3.90, 2.99 


3.24, 2.63 


237, 420 


12.5, 9.42 


12.77, 9.53 


1941, 1970 


0.3 


5 


2000 


1.795 


971 


14.0 


1988 


2.50, 2.00 


2.31, 1.86 


817, 950 


23.6, 16.1 


23.7, 16.1 


1986, 1990 


0.1 


3 


2000 


1.379 


496 


19.8 


1970 


1.782, 1.510 


1.573, 1.368 


433, 494 


32.1, 21.5 


31.6, 20.8 


1982, 1988 


0.03 


3 


2000 


1.192 


553 


32.9 


1969 


1.368, 1.234 


1.262, 1.163 


526, 555 


50.7, 33.2 


48.5, 31.1 


1991, 1994 


0.01 


2 


2000 


1.080 


286 


42.5 


1937 


1.179, 1.112 


1.097, 1.058 


273, 284 


65.6, 43.2 


58.4, 37.0 


1992, 1995 



Table II: Comparison of simulated and theoretical results. All values reported here are at or very close to steady state. 
Cells containing two numbers show theoretical predictions for p — 1/2 and 1, respectively. Domain size is L^ = Ly — 10"'' 
and I/z — 10^ (m). It is ensured that the domain is long enough so that the MTs don't reach the boundary for the given 
parameters. Here kgx ~ 5, Kst = 5, kgo = 0.5, Ksd ~ 500, KgT,gD = kgT,gDC. For nucleation rate we use Kn = Kgo, i.e., a 
D seed, except for the third line, where it is a T seed, Kn — KgT. jj^T/Nn is an estimate for the cap size - see Section fVI Bl 
Initially all Tu is free and its concentration is ctot- Notice that c < c^g. 



C. Non-steady-state phenomena 

Our numerical simulations are also capable of repro- 
ducing and explaining two other phenomena observed 
in experiments. Specifically, (i) in some cases the non- 
exponential MT length distribution and (ii) decaying 
oscillations in free tubulin concentration have been ob- 
served. 

MT length distributions. It is often mentioned in the 
literature that the steady state length distribution of 
MTs observed in the experiments is either exponential 
or bell-shaped H H ED, El E E^. Exponential dis- 
tribution agrees with our model. Inability to obtain a 
bell-shaped distribution seems to indicate a limitation 
to our model. While we don't exclude the possibility 
that some rates might depend on the MT length, as 
proposed by (a, 0, or on time spent in a given phase 
|44L l4^ . we suggest two alternative ways of obtaining 
bell-shaped distributions under certain conditions. First, 
one should be careful in determining when the system 
reaches the steady state in the experiment or simulation. 
As our simulations demonstrate, the system reaches the 
constant free tubulin concentration and MTs reach con- 
stant mean length relatively quickly. (The number of 
MTs does not change much henceforth.) By that time 
the MT length histogram is often bell-shaped as illus- 
trated in Fig. and Fig. This can be explained as 
follows. When MTs start growing from nucleation sites 
there is an excess of free tubulin. Therefore, the growth 
is originally unbounded leading to a Gaussian shape. If 



the cell edge (upper boundary) is far away, in the course 
of this growth the free tubulin concentration drops and 
reaches its steady state value. At this time the shape 
can still be close to a Gaussian (cf. (icL fig. 4]). This 
is followed by a long process of a shape change of the 
MT length histogram with free tubulin being constant. 
Eventually, this results in an exponential shape and in 
system reaching true steady state. This situation is well 
described in [4a,E3- 

A bell-shaped distribution can be also obtained in a 
bounded domain when the steady state concentration of 
free tubulin is high enough for unbounded growth of MTs 
if it were not for a cell edge (i.e., if c > c|^). In this situ- 
ation we predict positive exponential distribution of MT 
lengths, in the case when all MTs can reach identical 
maximal length restricted by the edge, consistent with 
[llL llfll|. As can be seen in experiments, however, MTs 
are curved and cell shape is not ideally spherical or circu- 
lar, so that different MTs experience different restrictions 
(e.g., 0,li3|)- This can lead to an MT length histogram 
of a bell-shaped form, in true steady state. 

Oscillations. In some of the simulations we observed an 
overshoot in free Tu concentration before the steady state 
was reached. Moreover, in some cases there were slight 
oscillations of free Tu concentration, as shown in Fig. 
Overshoots and large oscillations have been reported 
and modeled in the literature [S El IH H IH [Hi ii 
l52 |. It is believed that slow conversion of D- into T- 
tubulin in solution, after the depolymerization, is the key 
to understanding of such oscillations. These oscillations 
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Figure 5: Histograms of MT lengths after 80 (left) and 160 (right) seconds from the beginning of polymerization. The 
concentration of free Tu in both cases is close to steady state value: 7.63 and 7.51 ^M, respectively. Although the concentration 
has reached steady state after 80 sec, the length distribution of MTs is still changing. Here Kh = 10s~^, Ctot = W/J-M and rest 
of the parameters are specified in Table Hll 
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Figure 6: Histogram of MT lengths after 150 seconds from 
the beginning of polymerization. The concentration of free Tu 
is 1.382 I-lM and has reached its steady state value (to within 
random fiuctuations). Kh = 0.1s~ , ctot = SjiM and rest of 
the parameters are specified in Table El Note the bimodal 
character of the distribution. 



occur if initial free Tu concentration is sufficiently large. 
Free D-tubulin cannot polymerize. In our model we as- 
sume its conversion into T-tubulin to be instantaneous. 
We also use linear (and not higher order) dependence 
of nucleation rate on free Tu concentration, and a fixed 
number of nucleation sites. It is remarkable that under 
these restrictive assumptions the model produced some 
oscillations. We suggest the following explanation for 
their existence. If the hydrolysis is slow enough, the 
MTs grow quickly in the beginning, resulting in a large 



Figure 7: Small oscillations in the amount of polymerized 
Tu observed in numerical simulations for Kh = 0.1s~ , ctot = 
50/iM and A^„ = 10*. Other parameters are as in Table ITll 
Inset is a blowup of the region of oscillations. 



cap. When free Tu concentration changes quickly, the 
cap needs a relatively long time to adjust. This leads to 
a delayed response, and possibly to oscillations. Hence it 
might be that the ability to produce oscillations is inher- 
ent to MT structure, and that it can be magnified under 
certain experimental conditions. 



VII. CONCLUSIONS 

In this paper we analyze a novel model of MT dynamics 
in a domain bounded by the cell edge which involved 
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competition of individual MTs for tubulin. The model is 
based on a linear ID approximation of a MT structure. 
We consider the role of the boundary and extended the 
model to incorporate finite hydrolysis. 

One of the main results of the paper is an estab- 
lishment of a Hnk, by using analytical formulas, be- 
tween microscopic parameters describing polymeriza- 
tion/depolymerization and hydrolysis of individual units, 
and macroscopic (observable) characteristics of the MT 
dynamics and ensemble. We demonstrate how to ap- 
proximate macroscopic steady state behavior of MTs us- 
ing microscopic rates and vice versa, extract microscopic 
rates from macroscopic behavior. Hence, it is possi- 
ble to analytically and quantitatively predict the effect 
of changes in microscopic parameters on observable fea- 
tures, as well as to deduce microscopic changes from ob- 
served changes in macroscopic behavior, when relevant 
geometry and chemistry is taken into account. 

We consider the cell edge, and derive expressions char- 
acterizing MT distributions in the case of a finite pool 
of tubulin units. We perform Monte Carlo simulations 
to compare with our theoretical results, and find a good 
agreement between the two. 

The key ingredient in estabHshing a Hnk between 
micro- and macro-parameters is the cap model, which 
allows one to replace the actual cap consisting of many 
units by an effective single unit. The cap model behav- 
ior agrees with experiments measuring catastrophe fre- 
quency as a function of free tubulin concentration as well 
as with dilution experiments. 

We also use computational Monte Carlo model to pro- 
vide an explanation for a non-exponential MT length dis- 
tributions observed in experiments. 

We plan to develop in future a detailed model for 
an individual MT which will take into account the MT 
protofilament structure, and to use this model for investi- 
gating effects of various MAPs on the dynamic instability. 



each MT will be 1, which means that we find a maximal 
value among all K and then set At = 1/ iRaxMTs{K} 
[a technical note - we set At slightly below this value 
because the rate K for a given MT can change slightly 
as c is affected after each MT changes its length; these 
changes in c are usually very small]. Then, in general, 
each MT will have a chance to grow, shorten, or retain its 
length during At. We don't allow distribution of possi- 
ble number of length changes for an MT during At (only 
zero or one change is possible) . After the length of each 
MT is updated, we update c accordingly; after updating 
the lengths of all MTs, the hydrolysis cycle runs through 
all internal units of all MTs. The probability that a unit 
will hydrolyze during time At is taken as 1 — e^^'*'^*, 
assuming Poisson statistics. 

All MTs have their first unit in D state and this unit 
cannot be lost - this constitutes a simple nucleation seed 
with lower growth probability than when the MT tip is in 
T state (these units are not counted when calculating the 
lengths of the MTs) . Similarly, when the edge is relevant 
we can assume Ke = Kst- Such choices are made purely 
to reduce the number of parameters in the system, and 
are not essential for our purposes. 



Appendix B: EXPLICIT SOLUTION IN THE 
BOUNDED DOMAIN 



In what follows we obtain expressions for Nmt and m 
in the bounded domain at the steady state. Notice that 
in our model all the MTs of the maximum length L are 
technically in the growing phase, because their terminal 
unit can never become internal and therefore does not hy- 
drolyze. (These MTs cannot grow because of the edge.) 
The edge-induced catastrophe rate Ke will be governed 
by the smallest of the rates Kst, Kh- The discrete ver- 
sion of eqs. (I14|l and Ijl5|l determining a steady state at 
the boundary 
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Appendix A: COMPUTATIONAL MODEL 

In what follows we provide a short description of 
the numerical algorithm used for our simulations, which 
sHghtly differs from that used in [2j|- At time zero, MTs 
begin to grow from the nucleation seeds. At each simu- 
lation step, the time of this step is calculated by defining 
the rates of change in length (either growth or shorten- 
ing) for each MT. If the MT tip is in T (D) state then 
this rate of change \s K = Kgx(gD) + K^^xisD) ■ We de- 
mand that the maximal average number of changes for 



= K,Mg{L) - KgoMsiL - 1) - KsdMs{L - 1). (B2) 



After summing up these two equations we recover 

M,{L - 1) = {KfJ/K,D)Mg{L - 1), 



(B3) 



which is already known (cf. dJ) and so one of these 
equations is superfiuous. Another way to find Mg,s is to 
write the general solution of (tT4)l . l(T5)l : Mg = Ae'^^^ + 
B, M, - {KfJ /K,D)Ae-^/'' + {K°!^'/KgD)B and plug it 

in IHU yielding 5 = unless {K°^'K,D)/{KfJ Kgo) = 
1 in which case B is arbitrary. But this last condition 
implies A — > oo and hence Mg^g are constant inside the 
domain. At the lower domain boundary, eq. lf2ijl still 
holds. 
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The number of MTs is given now by 

L-l 



Nmt = J2(Mgil) + M,il)) + Mg{L) = 7V„ - TVo- (B4) 



1=1 



Using eqs. (QBl, C2|), GM ^nd I|B1|I and replacing the 
summation by the integration from to L we can deter- 
mine A and hence Nmt, which is given in eq. Ij28|l . 
Similarly, mean MT length is 



L-l 



m : 



Y,KMg{l) + Ms{l)) + LMg{L) 
1=1 

L-l 

Y^iMgil) + MS)) + M,{L) 



(B5) 



leading to eq. II29|I . 



1=1 
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